Disappearance of entanglement: a topological point of view 
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We give a topological classification of the evolution of entanglement, particularly the different 
ways the entanglement can disappear as a function of time. Four categories exhaust all possibil- 
ities given the initial quantum state is entangled and the final one is not. Exponential decay of 
entanglement, entanglement sudden death and sudden birth can all be understood and visualized 
in the associated geometrical picture - the polarization vector representation. The entanglement 
evolution categories of any model are determined by the topology of the state space and the dynam- 
ical subspace, the limiting state and the memory effect of the environment. Transitions between 
these types of behaviors as a function of physical parameters are also possible. These transitions are 
thus of topological nature. The symmetry of the system is also important, since it determines the 
dimension of the dynamical subspace. We illustrate the general concepts with a visualizable model 
for two qubits, and give results for extensions to A'^-qubit GHZ states and W states. 



PACS numbers: 03.65.Ud,03.65.Yz,03.67.Mn,02.40.Pc 



I. INTRODUCTION 

Quantum entanglement is widely accepted as a useful 
resource for communication and computation The 
preservation of this quantity through time is an im- 
portant goal in implementations of information transfer 
schemes and quantum computers. In most instances, en- 
tanglement is expected to decrease at long times because 
of the inevitable slow leakage of quantum information 
to the environment. However, recent work on models 
of decoherence of two entangled qubits has shown that 
the manner of this decay can be somewhat surprising. 
Besides the exponential decay of entanglement, i.e., the 
half-life (HL) behavior, it has been discovered that en- 
tanglement as a global property may abruptly terminate 
in a finite time, a phenomenon called entanglement sud- 
den death (ESD) [3. In subsequent work, ESD has been 
shown to be a rather general phenomenon: it can oc- 
cur when the environment is quantum, classical, Marko- 
vian, and non-Markovian [3-ilOi] . Oscillatory behavior of 
the entanglement function of time is observed in 

model calculations; this can take the form of entangle- 
ment sudden birth (ESB) if the two qubits are subject 
to a common bath [8|. Other non-monotonic evolutions 
of entanglement are also possible [III [13 • Very recently, 
the existence of ESD has been experimentally confirmed 
in optical and atomic systems 



13|. 



The elements of the density matrix are usually analytic 
functions of time, and the most typical behavior for them 
(or their envelopes) at long times is exponential decay. In 
ESD, in contrast, the entanglement measure goes to zero 
in a non-analytic fashion; this is because the typical en- 
tanglement measures are non-analytic functions of the 
elements of the density matrix. To date, we only have a 
"phenomenology" of possible behaviors of entanglement. 
The aim of this paper is to give a soundly based theoret- 
ical picture. In Sec. |TT1 we first categorize qualitatively 
the various possible time evolutions of the entanglement 
of two qubits and then show how the existence of these 



categories follows from the topology of the state space 
and of the spaces associated with dynamical evolutions, 
entanglement, and separability. 

For simplicity, we shall present formulas appropriate 
for the two-qubit case. However, the basic results gener- 
alize to N qubits. In Sec. IIIIl we illustrate the geomet- 
rical and topological arguments with a pure dephasing 
model for A^-qubit GHZ and W states. Finally, we sum- 
marize the results in Sec. IIVI 



II. TWO-QUBIT CASE 

A. state space and entanglement categories 

It is important to choose an appropriate representation 
of the state space. We use the polarization vector rep- 
resentation |14l4l6| . where the density matrix p, in the 
two-qubit case, is written as 
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i,j=I.X,Y,Z 



Here the components of the real 15-dimensional vector 
n are the expectation values of all physical observables, 
i.e., n = [{I ® X) , {I ®Y) , . . . {Y ® Z) , {Z ® Z)] where 
/ = (7/ is the 2x2 identity matrix, and X = ax , Y = , 
Z = az arc the Pauli matrices, n may be thought of as a 
generalization of the usual Bloch vector. This represen- 
tation is becoming increasingly popular to describe the 
results of experiments on multiple-qubit systems [13 ■ 
For iV qubits, the vector has dimM = 4''^ — 1 components 
(the number of elements of the su(2^) algebra). 

The state space M is the set of all physically admis- 
sible n, i.e., those that correspond to positive p. M is 
a compact, convex manifold with a boundary Bm- The 
fifteen matrices I ® X, etc., being generators of SU{4) 
group, form an orthonormal basis for the vector space in 
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which M is embedded and the inner product is chosen as 
(A, B) = Tr AS/4. The metric on M is the one that is 
induced by this inner product. 

Our measure of two-qubit entanglement is the Woot- 
ters' concurrence, denoted by C C is a continuous 
function on M that satisfies < C < 1: for separable 
states and 1 for maximally entangled states. The sepa- 
rable states defined by C{n) = form a set S that will 
play an important role below. is a real manifold with a 
boundary Bs- S C M and has "finite volume" in M, i.e., 
S is also 15-dimensional. Like M , S is compact and con- 
vex. Some points of the boundaries coincide: Bm n Bs 
is not empty and it forms a 4-dimensional manifold. S 
contains the origin n = 0, the totally mixed state. Im- 
portantly, it is known that S contains a ball of radius 
1 / centered at the origin j^O] . This is a lower bound 
for the radius rg of the maximum inscribed ball of sepa- 
rable states, and incidentally gives a connection between 
entanglement and purity. The chief difficulty in gener- 
alizing entanglement calculations to A^-qubit systems is 
that there is some degree of arbitrariness in all existing 
definitions of entanglement measure for N > 2. However, 
there is no arbitrariness in the definition of separability: 
for any N, a separable state is still any convex combi- 
nation of products of density matrices for the individual 
qubits, and any definition of entanglement must give zero 
on these states and no others. This is all that is required 
for our classification scheme. It is also still true for N 
qubits that M and S are compact and convex, and that 
dim M = dim S. The shape of M is complicated for large 
N but all pure states lie on its surface. The dimen- 
sionality of the submanifold of pure states is 2^+^ — 2 
and that of the pure separable states is 27V. 

The evolution of a quantum system is a smooth curve 
n(t) in M. This curve induces the continuous function 
C{t) = C{n{t)) that is the subject of this work. We shall 
take t S [0, oo) and only consider those trajectories with 
C(0) > and C(oo) = 0, i.e., those that begin in an 
entangled state and end in a separable state. (More gen- 
eral curves are certainly possible, and can also be usefully 
classified by the methods in this paper. For example, it 
is possible to give criteria for entanglement generation 
using our scheme.) We shall also assume the continuity 
as a function of time of all components of n{t), and all 
first derivatives. Let us denote the set of times when the 
entanglement is by Tq ee {t\C{t) = 0}. We define BSD 
(HL) behavior as any evolution such that Tq is of finite 
(zero) measure. 

Within these two larger classes we must also distin- 
guish subclasses because of the possibility of oscillations. 
We define four categories as shown in Fig. [T] Category 
A (approaching behavior) is defined by C{t) > Q for any 
finite t so that the curve of C (t) never actually hits the 
horizontal axis. Tq = in this case. This category in- 
cludes both monotonic and non-monotonic decay of en- 
tanglement [l^l- This is the "canonical" HL behavior. 
Category B (bouncing behavior) is defined by C{t) = 
only at isolated times for finite t so that Tq consists of 
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FIG. 1. Four categories of entanglement evolution. A: ap- 
proaching. B: bouncing. £: entering. O: oscillating. A and 
B belong to the HL class while £ and O belong to the BSD 
class. 

isolated points. Entanglement never quite dies for good 
in category B. Examples of entanglement evolutions in 
category B can be found in Tavis-Cummings model sys- 
tems when the initial state is in the "basin of attrac- 
tion" . Category £ (entering behavior) is defined by 
C {t) > ioT t <td < oo and C{t) = for t > td- This 
is the "canonical" ESD behavior. In category O (oscil- 
lating behavior) , Tq consists of disconnected intervals of 
finite width. Note these four possibilities are exhaustive 
given that the initial (final) state is entangled (separa- 
ble). The division into four categories is a result of two 
dichotomies: whether Tq is of finite measure and whether 
To is connected. Thus our classification is of topological 
nature. 

The existence of these categories can be explained by 
considering the trajectories n(t) and their relationships 
to S. Our basic premise is that three characteristics of 
any model determine the entanglement evolution cate- 
gories: the dynamical subspace D of the model, its lim- 
iting point Hoc = hmt_j.oo nit) and the memory effect of 
the environment. D is defined as the collection of possi- 
ble trajectories in a model. It is necessary to introduce 
D to understand how transitions between categories can 
happen. D can be of lower dimension than M if the evo- 
lution has any symmetries or if we limit the set of initial 
conditions in some way. The most important character- 
stic of any model is the set DCiS. If dim DOS ~ dim D, 
all four categories are possible. On the other hand, if 
dimD n S < dimD, only A- and i3-type behaviors are 
possible. This is because the measure of Tq in [0, oo) can 
be no greater than the measure oi D C] S in D. 

The limiting state noo exists for most physical decoher- 
ence processes; limit cycles and the like cannot be ruled 
out in general l2l| , but we will limit ourselves to the cases 
where S S and is unique. This constraint still leaves 
us with two possibilities: G Int(iS') (the interior of S) 
and noo G Bs- The first case guarantees the occurrence 
of categories £ and O while the second one could give rise 
to all four possibilities depending on other details of n{t), 
for example, whether noo is approached from or M\ S. 
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We take as a working definition that Markovian trajecto- 
ries satisfy the semigroup condition for all possible time 
partitioning [2^ . Then Markovian evolutions are ei- 
ther in categories A or £ while non-Markovian ones are 
typically in categories B and O. 



B. case study: a classical noise model 

The aforementioned three characteristics of a model 
do not uniquely determine the category of entanglement 
evolution. Transitions between categories are thus pos- 
sible by tuning some physical parameters of the model 
We now illustrate this with a general decoherence 
model that is nevertheless of low enough dimension that 
the topology can be visualized. 

The two-qubit Hamiltonian for the model is 

H = -\ \BKZK + bK{t)gK-SK]: 

K=A,B 

where Sk = {Xk,Yk, Zk)- The fields bA,B{t) are ran- 
dom time-dependent fields that decohere qubits A and 
B and they arc not correlated in time. Bk is a static 
field and is the noise coupling strength. For simplic- 
ity, we take Qa = Qb = 9- An important parameter in 
the model is 9, defined by cos6' = g ■ z. It is the angle 
between the energy axis and the noise axis. 9 = means 
that the noise is pure dephasing noise, while 9 = tt/2 
is transverse noise. If the correlation times of bA,B{t) 
are short (long) compared with l/gcos9, the system is 
typically Markovian (non-Markovian) [l^, HH . 

Decoherence in this classical noise model comes from 
the average over all noise histories bK(t). For more de- 
tailed discussions, see Ref. [1^1 and also the section on 
Bloch-Wangsness-Redfield theory in Ref. . 

We choose an initial state such that only nxx, nxY, 
nyx, nyy, and nzz ai'c nonzero, and this condition is 
preserved in the subsequent motion. We further require 
nxx = nyy and nxy = —nyx, which leaves only three 
independent parameters. This defines the dynamical sub- 
space D as a 3-dimensional slice of M. We refer to this as 
D3. D3 is large enough to accommodate essentially any 
decoherence dynamics given that (1) the two qubits are 
non-interacting; (2) noises on the two qubits are uncorre- 
latcd in time; (3) the effect of dephasing and relaxation 
can be separated; (4) the initial state is in D3. Thus 
D3 is the dynamical manifold of a rather general class 
of decoherence processes. Note that n-xx + "xy = 
is conserved in D3. The positivity of the density matrix 
requires 

02 > 2R^ + n\z < 3 (2) 
as > =^ nzz < 1 - 2i?^ and nzz > -1 (3) 
04 > 2R + nzz < 1 (4) 

where are the coefficients of the characteristic polyno- 
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FIG. 2. (Color online) Cross section of the dynamical sub- 
space D3 with nxY =0. a2 = is the ellipse, as = 
gives the parabola and the bottom of the isosceles triangle. 
04 = sets the two sides of the triangle. The entangled region 
is shaded where the filled arrows denote increasing direction 
of the concurrence. The green dashed line is the boundary 
of entangled and separable states, i.e., Bs n -D3. The fully 
mixed state is denoted by a red dot. This isosceles triangle 
corresponds to a projection of the tetrahedron in Ref. [2^. 

mial det(a;/ (g) I - p) = Ejlo(-l)^«r^^"^ 3- 

Applying these inequalities, we find that D3 is a cone 
in the (nxxinxyjnzz) coordinates. Its cross-section, as 
shown in Fig. [21 is an isosceles triangle with height 2 
and base 2 and the full manifold is generated by rotation 
about the nzz axis. The concurrence is given by 

C = max{0,R- {l + nzz)/^}. (5) 

Separable states S D form a set with a spindle shape 
on top and the entangled states form a torus-like shape 
on the bottom with triangular cross sections. The di- 
rection of the gradient of C is indicated by the filled ar- 
rows in Figl5] and the maximally entangled Bell states 
I*) = i (|01) +e*'^ |10)) with C = 1 form the circle 
R = 1 and nzz = — 1 on the lower surface of the cone. It 
is tempting to consider the minimum Euclidean distance 
to separable states in the polarization vector represen- 
tation as a geometric measure of entanglement [131 ■ 
It indeed works in but whether it qualifies as an en- 
tanglement measure in general is still an open question 

M- 

We thus fully characterized the entanglement topol- 
ogy of D3. Now we construct time evolutions in D3. 
The initial state is taken as the generalized Werner state 
w* = r I*) (^-l + (1 - r) 74/4 [a [l^- Without any loss 
of generality, we choose = in j^*), giving nxy = 0. 
This allows us to visualize the state and entanglement 
evolutions in a 2-dimensional picture. The initial state is 
"-(0) = (n-xx^nzz) = if,—r) and the state trajectory is 
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jiven by 



nxx {t) = {X®X) (t)=rC^^(i), 



-r e 



(6) 
(7) 



where = + Ff is the overah longitudinal re- 

laxation rate. = C'^C'^ describes overall dephas- 

ing process. Note C (0) = 1 and ( {oo) = if de- 
phasing occurs. The details of ({t) and Fi are model- 
dcpcndcnt. For Markovian noises, the Bloch-Wangsncss- 
Redfield theory applies and we have (t) = exp(— Ff^t), 
where Ff = Ff /2 -t- F^, Ff = g\sm^ OS'^ {Bt) /2 and 
F^ = g\cos^eS^{0)l2 ior K = A,B^. Here {uj) 
is the power spectrum function (Fourier transform of the 
noise autocorrelation function) of the classical noise pro- 
cess bK{t). For non-Markovian noises, non-cxponcntial 
behaviors in C,{t), such as damped oscillations, are possi- 
ble 0. 

The direction of the trajectory in is determined by 
the relative weight of dephasing and relaxation noise on 
the qubits. The trajectories are characterized by two pa- 
rameters, r which fixes the initial position, and 9, that 
weights the two types of noise. In the case of pure de- 
phasing [9 = 0), nzz is constant and the trajectory is 
horizontal with fioo = (0, — r). This is a visualization 
of decoherence free subspace where the constants of the 
decoherence dynamics can be used to encode informa- 
tion [sij . By contrast, when there is relaxation noise 
{9 7^ 0), the trajectory also moves vertically, so that 
when both types of noises are present, the trajectory will 
have = (0,0) € Int(5'). Only one trajectory with 
r = 1 and 9 = Q (pure dephasing), gives HL behav- 
ior since it remains in the entangled manifold except for 
noo = (0;~1) S Bs- All other trajectories show BSD. 
We may summarize the situation in topological language 
by saying that HL behavior requires that fioo € Bs- 

To pin down the precise category of the entanglement 
evolution, wc note that 



(8) 



*)/2. If the noise is Markovian 



C{t)^ 

where ^(t) = (1 — 
(non-Markovian) then C is monotonically decreasing (os- 
cillatory) d-i). In the latter case we find category O, as 
the state trajectory enters and leaves the separable re- 
gion as seen in Fig. [3^. In the case of Markovian noise, 
once the state leaves the entangled region, it leaves for- 
ever and wc find category £, as in Fig. ^p. There is a 
transition between this BSD behavior and HL behavior at 
the critical point r ~ 1 and 9 ~ Q. At this point wc have 
noo = (0,-1) S Bs- The transition is characterized by 
the fact that the "critical" trajectory intersects Bs n 
but not Int(S')n£'3. Markovian evolution yields category 
A behavior, as shown in Fig. [3]: while non-Markovian be- 
havior yields category B behavior as shown in Fig. [3ll. 
Another perspective is that the dynamical subspace for 
these trajectories with r = 1, = is the bottom disc D2 



of and dim_D2 n 5* = < dim_D2 = 2. The transition 
of To from the BSD class to HL class is thus a result of 
an topological transition of the dynamical subspace from 
A3 to D2. 



III. A/^-QUBIT CASE: GHZ AND W STATES 

Our examples have been drawn from the entanglement 
evolution of two-qubit systems. This is convenient, since 
the two-qubit concurrence is easily evaluated, and the 
lower dimensionality makes the examples relatively easy 
to visualize. However, it should be clear that the pre- 
cise definition of entanglement is not important for the 
topological categorization. Only the definition of sep- 
arability, which alone determines the set S, is crucial. 
As wc saw above, the topological properties of S, par- 
ticularly dim 5* = dimM, carry over to A^-qubit systems. 
There are indications that as the dimension increases, the 
percentage of the separable states in the physical states 
decreases f20l.[32l.[33|. 

The entanglement measure for general multi-partite 
quantum system is a topic of current research |34l436l | . 
The only easily computable entanglement measure for ar- 
bitrary A^-qubit mixed states is the negativity [36j , which 
measures the subsystem entanglement with respect to bi- 
partite partitioning of the whole system. It is defined as 



N{p) = 



\P 



1-1 



(9) 



Here the norm is taken to be trace norm and Ta denotes 
partial transpose on one of the two subsystems. N{p) cor- 
responds to the absolute value of the sum of the negative 
eigenvalues of p^^ , which according to Peres-Horodeck's 
criterion reveals the entanglement in p [37| . Here we use 
A^(/9) = as a working criteria for separability. 

We will use the negativity to illustrate the entangle- 
ment evolutions of A^-qubit GHZ and W states due to 
pure dephasing noise. The GHZ state and W states are 
of interest to the quantum computing community since 
they possess different types of multi-partite entanglement 
and are used in various protocols such as quantum secret 
sharing, tclcportation and super dense coding [35ll38l - l4l| . 
The A^-qubit GHZ and W states are defined to be 
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(10) 
+ 11...00)). 

(11) 



The dephasing process on a single qubit can be described 
in terms of two Kraus operators [12, 



|W) =_(|0...01)-l-|0...10) + 
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FIG. 3. (Color online) Concurrence evolution and state trajectory (inset) in the presence of Markovian and non-Markovian 
noise. The arrows denote time direction, (a) and (b) have both dephasing and relaxation present, while (c) and (d) are the pure 
dephasing case. The initial state is Bell state |^'+). We used ^'^^ = e"""^* cosO.Sf for non-Markovian noise and ^'^ 



for Markovian noise. For (a) and (b), T^ 



0.03. 



Note is smooth and characterize the dephasing pro- 
cess in the xy plane in the Bloch vector representation, 
just as in D^. For simplicity, we assume that all qubits 
are independent and they share the same dephasing func- 
tion C{t). Thus the full dynamics of the A^-qubit system 
can be described by 



Pit) 



E 



E,,E,,---E,,poEI---EIeI. (13) 



For the negativity, we make the {1){N — 1) bipartite 
partition, i.e., the partial transpose is applied to one of 
the N qubits. Due to the permutation symmetry of GHZ 
and W states, as well as the iV-qubit Kraus operators, it 
docs not matter which qubit is picked, indicating that the 
entanglement is evenly distributed among the TV qubits. 
Given these assumptions, the negativity of the iV-qubit 
GHZ and W states are given by 



A^GHZ =|Cl'W/2 



w 



N 



(14) 
(15) 



Notice that due to the smoothness of entanglement 
evolutions in categories £ and O cannot happen. This 
can be explained by the topological argument as follows. 

For both the GHZ and W states, the dephasing process 
does not expand the Hilbert space. For the A^-qubit GHZ 
states, the density matrix can be expanded with identity 



and elements of the su{2) algebra 



Pghz 



i=X.Y,Z 



(16) 



where ai are defined on the two cat states, for example, 
ax = |0)®^(1|®^ + |1)®^(0|®^. In other words, the 
dynamical subspace -Dghz for GHZ states under dephas- 
ing is a 3-dimensional Bloch ball, as seen in Fig. |4l The 
negativity is given by 



^(pGHz) 



(17) 



The separable states in terms of negativity live on the 
vertical symmetry axis. We have 

dim -Dghz H 5 = 1 < dim -Dghz = 3 

thus only entanglement evolutions in categories A and B 
are possible. 

For the W states, the density matrix is expanded with 
identity and elements of the su{N) algebra. Thus the dy- 
namical subspace -Dw is N'^ — 1 dimensional. However, 
when calculating negativity, the partial transpose intro- 
duces N extra bases into the Hilbert space and p^^{t) 
contains elements in su(2N) algebra. The partial tran- 
pose moves N — 1 elements f3i and their complex con- 
jugates /3* in the original density matrix to the blocks 
of the newly introduced bases. The positions of these 



6 





L — 


N =o\ 

"ghz " \ 


















t 1 
N =1/2/ 

GHZ / 



FIG. 4. (Color online) Effective Bloch sphere representation 
for the dynmaical subspace Dghz ■ The north and south 
pole are either |0)^ or States resting on the red line 

connecting north and south pole have zero negativity. States 
on the green equator are maximally entangled. 



elements /3i are related to the partial transpose: sup- 
pose the partial transpose acts on the fc'th qubit, then fii 
are the coefficients of |. . . . . .) (. . . 1 . . .|, where and 1 
label the fc'th qubit. Take the 3-qubit case for example, 
the component /3i|001)(100| is mapped into /3i|101)(000|, 
where the partial transpose is taken on the first (left- 
most) qubit. On the other hand, components such as 
|001)(010| or |100)(100| are left unchanged under the ac- 
tion of partial transpose on the first qubit. 

The negativity for any density matrix expandable by 
the W state bases is given by 



1 

TV' 



(18) 



Note that the condition of separability N{pw) = elim- 
inates 2{N — 1) degrees of freedom. Thus 



dim(i:>w n 5) = (TV - if < dimL> 



w 



1 



and we conclude that entanglement evolutions in cate- 
gories £ and O are not possible. 



IV. SUMMARY 

In summary, we found that there are four and only 
four types of natural behaviors for the time evolution of 
entanglement, given the initial state is entangled and the 
final one separable. These categories are determined by 
the dimensionality and intersection properties of sets in 
M. Since these properties are preserved by continuous 
deformations, they are topological. Three characteristics 
determine the categories: the dynamical subspace D of 
the model, the limiting state noo and the memory effect 
of the environment. Category A occurs (typically) for 
Markovian systems when rioo S Bs, category B occurs 
(typically) for non-Markovian systems when rioo G Bs- 
Category £ occurs (typically) for Markovian systems, 
while category O occurs (typically) for non-Markovian 
systems. For those two categories, noo can be either on 
the boundary or in the interior of S. 

Model studies of entanglement evolution have shown a 
wide variety of behaviors, and a unifying picture of these 
behaviors has been lacking. The topological approach 
given here provides such a picture. The qualitative be- 
havior is determined by relative dimensions of the dy- 
namical subspace, the space of separable states, and the 
space formed by their intersection. The determination 
of these dimensions, together with the location of the 
asymptotic point, tells us when transitions between dif- 
ferent types of evolution can be expected. Precise means 
of computing the dimensions, including the important 
role of symmetry, will be given in a future publication. 
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